Exactly-solvable coupled-channel potential models of atom-atom magnetic Feshbach 
resonances from supersymmetric quantum mechanics 
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Starting from a system of N radial Schrodinger equations with a vanishing potential and finite 
threshold differences between the channels, a coupled N x N exactly-solvable potential model is 
obtained with the help of a single non-conservative supersymmetric transformation. The obtained 
potential matrix, which subsumes a result obtained in the literature, has a compact analytical form, 
as well as its Jost matrix. It depends on N(N + l)/2 unconstrained parameters and on one upper- 
bounded parameter, the factorization energy. A detailed study of the model is done for the 2x2 case: 
a geometrical analysis of the zeros of the Jost-matrix determinant shows that the model has 0, 1 or 2 
bound states, and or 1 resonance; the potential parameters are explicitly expressed in terms of its 
bound-state energies, of its resonance energy and width, or of the open-channel scattering length, 
which solves schematic inverse problems. As a first physical application, exactly-solvable 2x2 
atom-atom interaction potentials are constructed, for cases where a magnetic Feshbach resonance 
interplays with a bound or virtual state close to threshold, which results in a large background 
scattering length. 
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I. INTRODUCTION 



Coupled-channel quantum-scattering models are ex- 
periencing today a strong renewal of interest, mainly 
thanks to the impressive experimental progress in the 
field of ultracold gases. Actually, an indispensable tool 
to control the atom-atom interactions in these systems 
relies on the coupling between channels defined by dif- 
ferent hyperfine states of the atom-atom pairs. When 
dipped in a magnetic field, these hyperfine states have 
threshold energies which vary linearly with the field. 
At ultracold temperatures, the lowest-threshold state is 
the only open channel, while states with higher thresh- 
olds are closed. When varying the field in a well- 
chosen range, quasi bound states of the closed chan- 
nels can appear as resonances in the open channel, a 
phenomenon known as a "Feshbach resonance" , a phe- 
nomenon first studied in nuclear physics [I], When 
a Feshbach resonance crosses the open-channel thresh- 
old, due to magnetic-field variation, the scattering length 
a, which effectively controls the atom-atom interaction, 
goes through infinite values, switching from positive to 
negative sign [!, 0, [f| . This very spectacular phenomenon 
is now known as a magnetic-field-induced Feshbach res- 
onance, or just "magnetic Feshbach resonance" . 

The practical importance of magnetic Feshbach reso- 
nances has motivated various theoretical models, which 
can be classified in three categories: (i) microscopic 
models, which should in principle deduce magnetic- 
Feshbach-resonance properties from many-electron cal- 
culations; (ii) effective potential models, which reduce 
the complexity of the many-electron problem to a two- 
atom coupled-channel problem (usually two channels arc 
enough), where the interaction between the two atoms is 
modeled by a symmetric (two by two for a two-channel 
model) potential matrix; (iii) effective scattering-matrix 



models, which reduce the role of the underlying interac- 
tions to its impact on the atom-atom open-channel scat- 
tering matrix. Present-day theoretical descriptions of ul- 
tracold gases only require the knowledge of the atom- 
atom scattering length, which is directly related to the 
open-channel scattering matrix. As far as practical ap- 
plications are concerned, the third category of models is 
thus sufficient. There, a magnetic Feshbach resonance 
is described as a pole of the scattering matrix in the 
complex wave- number planes, like any resonance [||, and 
the whole complexity of the many-electron or atom-atom 
problem is reduced to a few parameters of the scattering- 
matrix Pade expansion Q . These parameters can be nu- 
merically obtained, e.g., with the help of the reaction- 
matrix method [8| , from a given microscopic or effective- 
potential model. 

In several contexts, the use of effective-scattering- 
matrix models is however progressively felt to be insuffi- 
cient. This is for instance the case when a magnetic Fesh- 
bach resonance occurs with a large background scattering 
length, due to a bound or virtual state in the open chan- 
nel, close to its threshold Q ■ Such an open-channel state 
is also called a "potential resonance" because it natu- 
rally occurs in a potential model, even in a single-channel 
case. Other situations where a more detailed knowledge 
of the atom-atom interaction than just the scattering 
length might be necessary are cases where molecules can 
be formed, as in crossovers between a Bardeen-Cooper- 
Schieffer supcrfluid and a Bose-Einstein condensate, or in 
Bose-Einstein-condensate collapses. None of these cases 
probably requires to solve the full many-body electronic 
problem; effective-potential models, on the other hand, 
look like a reasonable approximation, as they allow for 
a realistic description of the atom-atom interaction in 
terms of the accessible channels and as a function of the 
radial coordinate r between atoms. 

There is thus an interest for exactly-solvable coupled- 
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channel potential models with threshold differences. The 
first example coming to mind is probably the coupled 
square-well potential, which can display both potential 
and Feshbach resonances, as well as bound states @. 
This model has however two drawbacks: first, despite its 
simplicity and exactly-solvable character, its scattering- 
matrix poles are given by rather complicated implicit 
equations. Second, its discontinuous form factor is rather 
limitative and very different from the known long-range 
atom-atom polarization interaction. The next choice to- 
wards realistic atom-atom interactions is thus a purely 
numerical resolution of the coupled- channel Schrodinger 
equation with smooth phenomenological potentials. This 
lack of exactly-solvable potentials can be related to the 
poor knowledge of the scattering inverse problem (i.e., 
the construction of a potential in terms of its bound- 
or scattering-state physical properties) in the coupled- 
channel case with threshold differences [10| . In Ref. [TO ] 
however, an exactly-solvable coupled-channel potential 
with threshold differences is derived, two remarkable fea- 
tures of which are the compact expressions provided both 
for the potential and for its Jost matrix. Since the Jost- 
matrix completely defines the bound- and scattering- 
state properties of a potential model fl2l , [T^ . such an 
analytical expression seems very promising in the con- 
text of the scattering inverse problem. 

The work of Cox has however received little attention, 
probably because it is plagued by two problems. First, 
the way of getting the potential is rather complicated and 
mysterious: the paper mostly consists in a check that the 
provided analytical expression for the solutions satisfies 
the coupled-channel Schrodinger equation with the pro- 
vided analytical expression for the potential. Not much 
information is given on how these expressions were ob- 
tained, which makes any generalization of the method 
impossible. The second problem, already stressed in 
Ref. [ll|, is that, despite the compact expression of the 
Jost matrix, calculating the corresponding bound- and 
resonant-state properties is a difficult task because these 
states correspond to zeros of the determinant of the Jost 
matrix in the intricate structure of the energy Riemann 
sheet, which has a multiplicity 2 N for N channels. 

The first problem was solved recently, when it was re- 
alized that the Cox potential, at least in its simplest form 
(q = 1 in Ref. fllj), can be obtained by a single super- 
symmetric transformation of the zero potential [lil . Il5| . 
This leads to a much simpler derivation of this poten- 
tial and naturally enables several generalizations of it; 
in particular, the initial potential is now arbitrary. The 
transformation used to get this result belongs to a cate- 
gory of supersymmetric transformations not much used 
up to now, namely transformations that do not respect 
the boundary condition at the origin (the so-called non- 
conservative transformations, see e.g. Ref. [Til]): a so- 
lution of the initial potential vanishing at the origin is 
transformed into a solution of the transformed poten- 
tial which is finite at the origin. This feature makes the 
transformation of the Jost matrix more complicated to 



calculate than for usual conservative transformations but 
it is also the key to get potentials with non trivial cou- 
pling. In Sec. |TT] below, we give several alternative ex- 
pressions for the Cox potential and explicitly make the 
link between the supersymmetric derivation and the ex- 
pressions found in Ref. II], for an arbitrary number of 
channels N. With respect to Refs. [3, EH, this result 
is new as these references mostly concentrate on general- 
izations of the Cox potential allowed by supersymmetric 
quantum mechanics and on N = 2 examples. We also 
show in Sec. [II] that the Cox potential contains the max- 
imal number of arbitrary parameters allowed by a single 
non-conservative supersymmetric transformation, which 
makes it the most interesting potential from the point of 
view of the scattering inverse problem; we also derive a 
new necessary and sufficient condition for the regularity 
of the potential. 

The second problem, i.e., the calculation of bound- and 
scattering-state properties from the analytical Jost func- 
tion, is touched upon in Sec. IIIIl There, the discussion 
is limited to N = 2 (a case complicated enough from 
the mathematical point of view but very rich already 
from the physical point of view). First, the number of 
bound states and resonant states is studied geometri- 
cally in terms of the potential parameters, as well as 
the necessary and sufficient condition for a regular po- 
tential (several mistakes made in Ref. in particular 
regarding the number of bound states, are corrected in 
passing). Second, schematic inverse problems are solved, 
where the potential parameters are expressed in terms of 
physical quantities like bound-state energies, resonance 
energy and width; third, the low-energy behavior of the 
open-channel scattering matrix is studied, with ultracold 
gases in mind. This discussion makes possible a first 
practical use of this potential as a schematic model of 
atom-atom magnetic Feshbach resonances, which is de- 
scribed in Sec. IIVI There, an exactly-solvable model is 
established in cases where a magnetic Feshbach resonance 
interplays with a potential resonance, which results in a 
large background scattering length, either positive (in- 
terplay with a bound state) or negative (interplay with a 
virtual state). This physical context is mostly inspired by 
Ref. 0]. Sec. fVl finally summarizes our findings and dis- 
cusses possible extensions of them, in particular to other 
fields of physics where coupled-channel models are known 
to play an important role. 



II. THE COX POTENTIAL FROM 
SUPERSYMMETRIC QUANTUM MECHANICS 

Let us first summarize the notations used below for 
coupled-channel scattering theory 0, [l2|, EH- We con- 
sider a multichannel radial Schrodinger equation that 
reads in reduced units 

Hip(k,r) = K 2 iP(k,r), (1) 
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with 



H 



dr 2 



V, 



(2) 



where r is the radial coordinate, V is an N x N real 
symmetric matrix, and ip may be either a matrix-valued 
or a vector-valued solution. By k we denote a point in 
the space C^, k = {ki, . . . , fcjv}, ki S C. A diagonal 
matrix with non-vanishing entries /c^ is written as if = 
diag(/c) = diag(fci, . . . , kjy)- The complex wave numbers 
ki are related to the center-of-mass energy E and the 
channel thresholds A 1; . . . , A^v, which are supposed to 
be different from each other, by 



E 



(3) 



For simplicity, we assume here that the different channels 
have equal reduced masses, a case to which the general 
situation can always be formally reduced [Hj]. We also 
assume potential V to be short-ranged at infinity and to 
support a finite number M of bound states. Under such 
assumptions, the Schrodinger equation has two N x N 
matrix- valued Jost solutions which allow one to construct 
the Jost matrix F{k) defining both scattering and bound- 
state properties. The scattering matrix, which is sym- 
metric, reads 



S(k) = K- l l 2 F{-k)F- l {k)K 1 ' 2 

= K^ 2 [F- 1 (k)} T F T (~k)K- 1 ^ 2 , 



(4) 



with T meaning transposition and —k = 
{— ki km}. The zeros of the determinant of 
the Jost matrix, which are defined by detF(fc) = 0, thus 
correspond to poles of all the elements of the scattering 
matrix. Bound states correspond to such zeros k m , with 
m = 1, . . . , M, lying on the positive imaginary fcj axes for 
all channels: k m i — in m i with K m i > and i = 1, . . . , N. 
The corresponding energies, E m = -k^ + A.;, lie below 
all thresholds. For simplicity, we call virtual state any 
other zero of the Jost-matrix determinant corresponding 
to a real energy below all thresholds, but not lying on all 
the positive imaginary ki axes. Finally, we call resonance 
any zero of the Jost-matrix determinant not lying on 
the imaginary ki axes, hence corresponding either to a 
complex energy or to a real energy above at least one 
threshold. Note that for a resonance to have a visible 
impact on the physical scattering matrix it should be 
located sufficiently close to the real axis. 

Let us then summarize the main results from super- 
symmetric quantum mechanics in the coupled-channel 
case [16|, llli- Starting from an initial potential V and 
its solutions ip, a supersymmetric transformation allows 
the construction of a new potential 



V(r) = V{r) - 2U'(r) 



with solutions 

i>(k,r) = 



- T + U{r) 
dr 



ip(k,r) 



(5) 



(6) 



where the so-called superpotential U is expressed 
terms of a square matrix a by 

U(r) = a' '(r)a- 1 (r) . 



m 



(7) 



Matrix a is called the factorization solution; it is a solu- 
tion of the initial Schrodinger equation 



Ha(r) = -)C 2 a(r) , 



(8) 



where K, — diag(rt) = diag(K 1; . . . ,/cjv) is a diagonal ma- 
trix called the factorization wave number, which corre- 
sponds to an energy £ lying below all thresholds, called 
the factorization energy. The entries of K, thus satisfy 



S 



Aj; by convention, we choose them positive: 



Ki > 0. Equation ^ implies that all the physical prop- 
erties of the transformed potential can be expressed in 
terms of those of the initial potential, in particular its 
Jost matrix and scattering matrix. 

Let us now apply these results to a vanishing initial 
potential 1^ = 0, for which the Jost matrix and scatter- 
ing matrix are identity, S(k) = F(k) = I. For a given 
factorization energy, the most general real symmetric su- 
perpotential depends on an TV-dimensional real symmet- 
ric matrix of arbitrary parameters, i.e., on N(N + l)/2 
real arbitrary parameters [l5|]. When V = 0, the corre- 
sponding factorization solution can be written as 

cr(r) = cosh(«r) + JC^ 1 smh(nr)Uo (9a) 
= (2/C)- 1 [exp( K r)(/C + U ) + exp(- K r)(/C - U Q )], (9b) 

which ensures that the resulting potential V is regular at 
the origin, and where the arbitrary parameters explicitly 
appear as the value of the (symmetric) superpotential at 
the origin, Uq = U(0); exp(±Kr), cosher) and sinh(Kr) 
are diagonal matrices with entries exp(±Kir), cosher) 
and sinh(Kir) respectively. According to Ref. [l5l ]. when 
K, + Uq is invertible, the transformed Jost matrix reads 



F(k) = {K- iK)- 1 ^ - iK). 



(10) 



This is the Jost function obtained by other means in 
Ref. ll| in the case q = 1. However, it was not re- 



alized there that the corresponding potential could be 
simply expressed in terms of a solution matrix a, using 
Eqs. (O and (0). In that reference, a compact expression 
for the potential is found [see Eq. (JT5J) below] but writ- 
ing §§§ and ([7]) is much more elegant because both the 
potential ([5]) and its Jost function (fT0|) are expressed in 
terms of the same parameter matrix Uq. Nevertheless, 
this procedure also presents several disadvantages: cal- 
culating the potential requires several matrix operations 
(inversion, product, derivations); moreover, the parame- 
ters in Uq should be chosen so that the factorization solu- 
tion is invertible for all r, a condition not easily checked 
on Eqs. ©. 

Let us now derive an alternative form for the factoriza- 
tion solution, which solves both these inconveniences. In 
Ref. [TH, the possibility of rank (tC + U ) < N in Eq. ([90) 
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has been studied, which leads to an interesting asymp- 
totic behavior of the superpotential but which reduces 
the number of parameters in the model. Here, in order 
to keep the maximal number of arbitrary parameters in 
the potential, we choose K. + Uq invertible. The factor- 
ization solution (I9bjl can then be multiplied on the right 
by 2(JC + Uo)^ 1 ^ 1 / 2 , which leads to the factorization so- 
lution 

a{r) = /CT 1/2 [exp(nr) + exp(-«r)Xo] . (11) 

According to Eq. (J7J), the superpotential, and hence the 
transformed potential, is unaffected by this multiplica- 
tion. The symmetric matrix Xq now contains all the 
arbitrary parameters. The link between the two sets of 
parameters is given by 

X Q = /C- 1/2 (/C-C/ )(/C + C/ )" 1 /C 1/2 , (12) 
U = JC^il-Xo^I + Xo)- 1 ^ 2 . (13) 



The last expression is exactly equivalent to Eq. (4.7) of 
Ref. [IH for q = 1, which reads 

V(r) = 2e~ Kr [I - A(2K.)- 1 e- 2,ir y 1 (AK + K.A) 

x [/ - e - 2Kr (2IC)- 1 Ay 1 e- Kr , (18) 

provided one defines matrix A as 

A = -2K>l 2 XoK>l 2 

= -2(IC-Uo)(IC + Uo)- 1 IC. (19) 

The second advantage of writing is that it eas- 
ily leads to a necessary and sufficient condition on the 
parameters to get a potential without singularity at fi- 
nite distances. This condition is positive definiteness of 
matrix I + Xq : 

I + X Q >0. (20) 

The potential has a singularity when a(r) is non- 
invertible, i.e., when det[/ + A(r)] vanishes for some 
r. Using Eq. ([To]), we find that this is equivalent to 
the existence of tq > such that dety(ro) = with 
Y(r) = exp(2rer) + Xq. Assume now that detF(r) ^ 
Vr > 0. Since detF(r) = J^[ i=1 y%{r) where yi(r) are the 
eigenvalues of Y(r), we conclude that yi(r) ^ for all 
i = 1, . . . , N and r > 0. But since for sufficiently large r, 
Xq becomes a small perturbation to exp(2Kr), all eigen- 
values of Y(r) should be positive for r > and in par- 
ticular at r — 0, thus proving the necessary character of 
the above condition. 



Equation (jlip can also be written as 

o-(r) =/CT 1/2 [/ + A(r)]exp(Kr), (14) 

where 

X(r) = exp(— kt)Xq exp(— nr). (15) 

With respect to writing (|5a|) and (|9bp . Eq. Q14j) presents 
several advantages. First, it allows for a simple calcula- 
tion of the superpotential 

U(r) = /C-2/C 1/2 X(r)[/ + X(r)]- 1 /C 1/2 

= -/C + 2/C 1/2 [7 + X(r)]- 1 JC 1/2 . (16) 

The last expression is particularly convenient since the r 
dependence is limited to one factor of the second term; 
the potential can thus be explicitly written as 



(17) 

I 

The sufficiency follows from the observation that Y(r) 
is positive definite for any r > 0, together with Y(0) = 
I + Xq. Indeed, if Y(r) is positive definite, the in- 
equality (q\Y(r)\q) > holds for any q £ Zjy Here 
(p\q) = J^iLiPiQi i s the usual inner product in the N- 
dimensional complex linear space Ljv, with pi, qi being 
coordinates of the vectors p,q G Ln with respect to an 
orthonormal basis. But since (q\Y(r)\q) = (q\X \q) + 
(q\exp{2Kr)\q) > (q\X \q) + (q\q) = (q\X + I\q) [we re- 
call that r > 0, Ki > and exp(«r) is a diagonal matrix 
with entries exp(/«jr)], positive definiteness of I + Xq im- 
plies positive definiteness of Y(r) for r > 0. 

Having established this condition on Xq, one can get 
the condition in terms of Uq, using Eq. (|12[) . Since 

I + Xq = 2K 1/2 {JC + Uq)- 1 K 1/2 1 (21) 

the necessary and sufficient condition to get a regular 
potential is positive definiteness of matrix }C + Uq\ 

K + U >0. (22) 

Since the (diagonal) elements of K, are positive and in- 
crease when the factorization energy decreases, this con- 
dition has a simple interpretation: it just puts some up- 
per limit on the factorization energy. 

Finally, Eq. p^l) shows that the condition det A ^ 
required in Ref. [ill] is not required here. In Cox' paper, 
this condition does not appear in the potential expres- 
sion, which is valid in the general case, but only in the 
derivation of the proof; the fact that this condition is not 



I 

V(r) = 4/C 1/2 [/ + A(r)]- 1 A w (r)[/ + A(r)]- 1 /C 1/2 

= -4/C 1/2 (e Kr + Xoe-^y 1 {XqK. + KXq) (e Kr + c'^Xq) ~ l K 1/2 . 
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required here illustrates the efficiency of the supersym- 
metric formalism. Equation (|19[) also implies that rank 
(fC + Uo) < N corresponds to det A = oo, a case also not 
considered in Ref. The supersymmetric treatment, 
on the contrary, allows this case [l5| ; our approach 
thus subsumes the results of Ref. [llj in several respects. 



III. GENERAL PROPERTIES OF THE 2x2 COX 
POTENTIAL 

Having established a connection between the Cox po- 
tential and supersymmetric quantum mechanics, we now 
proceed in a more detailed analysis of its properties for 
the simplest particular case N — 2. As is happens, this 
case is not only complicated enough to deserve a ded- 
icated analysis, but also sophisticated enough to make 
the solution of several interesting inverse problems pos- 
sible. 



A. Explicit expression of the potential 

For N = 2, the arbitrary parameters entering the Cox 
potential are the entries of the superpotential matrix at 
the origin, 

/ oi P \ 
U ee 1/(0) = , (23) 

\ p a 2 J 

I 



and the factorization energy £. The corresponding fac- 
torization wave number, n = (ki,k 2 ), is made of two 
positive parameters K\ and K2 which are not indepen- 
dent of each other: they should satisfy the "threshold 
condition" [see Eq. ©] 

k\ - k\ = A. (24) 

Here and in what follows we put for convenience Ai = 0, 
A 2 = A > 0. 

In terms of these parameters, the necessary and suffi- 
cient condition for a regular potential, i.e., JC+Uq positive 
definite, can be written for instance 

Ki > — oti, (25a) 
2 

K2 > — a 2 . (25b) 

Ki + ai 

This puts an upper limit on the factorization energy in 
terms of the parameters appearing in Uq [see Eq. (|43|) 
and Fig. [5] below] . 

Two explicit expressions for the superpotential are 
given in Ref. (lif . Using Eqs. ([5]) and ([To]) , one gets what 
is probably the simplest possible explicit expression for 
the potential itself: 



2Kir xnKi + [2X11X22K1 - x\ 2 (n x + k 2 )] e 2fi2r + x 2 2 (xnX22 - xi 2 ) K ie 



[1 



xue- 2 ^ r + x 2 2e~ 2K2r + (xnx 22 - x\ 2 ) e -^i+^)r 



U12 = -4xi 2 V«i«2e 



— (ki+K2)i- 



ki + k 2 + xu(k 2 - Ki)e" 2Kir + x 22 ( Kl - K 2 )e~ 2K2r - (x n x 22 - x\ 2 ) («n + K 2 )e~ 2< - Kl+K ^ r 
[l + xne- 2 ^ + x 2 2e- 2K2r + (xnx 2 2 ~ x\ 2 ) e -^i+^)r^ 2 



(26a) 



(26b) 



The element £22 is obtained from Eq. (|26a[) by the re- 
placement ki <-> K2 and xn ^> X22- Here, we have used 
the symmetric matrix 



Then, according to Eq. ([10]) , the Jost matrix for the Cox 
potential reads (see also Refs. [Ill HiL Hll]) 



Xll X12 



X12 X22 



(27) 



which is related to matrix J23J) by Eqs. H2J) and (|13|). 
In the following, as we are mostly interested in the Jost- 
matrix properties, we shall rather use matrix Uq. 
B. Zeros of the Jost-matrix determinant 



F(k,p) = 



03 



k-\-iai 

p+i«2 p-\~tK2 



(29) 



Let us denote for convenience the channel wave num- 
bers as k\ — k and k 2 — p, with the threshold condition 



k 2 -p 2 = A. 



(28) 



The determinant of the Jost matrix coincides with 
the Fredholm determinant of the corresponding integral 
equation (T^ |; it reads here 
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Figure 1: Geometrical representation of Eqs. (|34ap (first column, solid lines) and (|34b[) (first column, dashed lines), 
and positions of the corresponding roots of system (|3ip in the complex k (second column) and p (third column) planes. 
Various values of the parameters oti, 012 are chosen, which imply various numbers of bound, virtual and resonant states: 

(a) ot\ < 0, 012 < — VA, two bound states (star and diamond), two virtual states (square and triangle), no resonance; 

(b) ai > 0, a.2 < — vA, one bound state (star), one virtual state (square), appearance of a resonance (diamond); 

(c) Qi > 0, Q2 > 0, no bound state, two virtual states (star and square), one resonance (triangle and diamond); 

(d) Qi > 0, 0:2 > \/A, no bound state, no resonance, four virtual states. Increase of either ai or a.2 leads to: (a) (thin 
dashed lines) disappearance of a bound state; (b) appearance of the resonance. 



the system of equations for finding the zeros of f(k,p), 



f(k,p) = detF(k,p) 



{k + ia 1 ){p + ia 2 ) + /3 2 
(k + ini)(p + 



(30) 



fc 2 -p 2 = A, (31a) 
{k + i ai )(p + ia 2 ) +f3 2 = 0, (31b) 



The zeros of this Jost determinant in the k and p com- 
plex planes, which correspond to bound, virtual or reso- 
nant states, are functions of the parameters a\, a 2 and f3 
only. From here and the threshold condition (|28p , follows 



which is equivalent to the fourth-order algebraic equa- 
tion: 



k' + ia±k + a 2 k + ia^k + 04 = , 



(32) 
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where 



ai 


= 2ai, 






(33a) 


02 


2 2 
= a 2 - a x 


-A, 




(33b) 


as 


= 2[ ai (a 2 2 


-A) 


-a 2 /3 2 ], 


(33c) 


0.4 


2/ 2 

= -Oi(a 2 


-A) 


+ 2a 2 p 2 a 1 -p 4 . 


(33d) 



We notice that after substitution k = iX, Eq. (f3"2"| be- 
comes an algebraic equation in A with real coefficients. 

Its four roots are thus either real numbers, which corre- 
spond to real negative energies (bound or virtual states), 
or mutually-conjugated complex numbers, which corre- 
spond to mutually-conjugated complex energies (reso- 
nant states). Basing on this property, we will use in 
what follows a geometric representation of the system of 
equations which allows for a visualization of the zeros of 
f(k,p) in the parameter space. 

Let us first consider bound and virtual states, which 
correspond to solutions of system pip with k and p 
purely imaginary. After substitution k = iX, p = ip, 
with A and p real, these equations define two hyperbolas 
in the (A, p)-plane, 



p 2 -A 2 = A, 

(X + a 1 )(p + a 2 ) =/3 2 , 



(34a) 
(34b) 



the positions of which are defined by the values of the 
parameters ai, a 2 , /3 and A. 

The roots of system (13"4"|) that correspond to bound 
and virtual states are the intersection points of these hy- 
perbolas. Different possibilities of hyperbola locations 
are shown in Fig. [1] The solid-line hyperbola corre- 
sponds to the threshold condition (|34a|) ; its semi-major 
axis is y/~K and its slant asymptotes are given by p = ±A. 
The dashed-line hyperbola corresponds to Eq. ()34b|) ; its 
asymptotes are given by A = — a± and p = —ct2- The 
abscissa (resp., ordinate) of a crossing point in the (A, p)- 
plane gives the position of the corresponding zero on the 
imaginary axis in the fc-plane (resp., p-plane), as shown in 
the second (resp., third) column of Fig. Q] Bound states 
correspond to X, p > 0, i.e., to intersection points lay- 
ing in the first quadrant of the (A, p)-plane, while virtual 
states correspond to intersections in the second, third 
and fourth quadrants. It is clearly seen on Fig. [1] that 
the two hyperbolas (|34ap and (|34b[) cross in cither two 
or four points. Moreover, they can have zero, one or 
two intersections in the first quadrant, which means that 
the potential has either zero, one or two bound states. 
This contradicts Ref. [HI], where it is said that the po- 
tential does never support bound states. Since Eq. (|3^|) 
is fourth order, when the hyperbolas cross in four points, 
the Jost determinant does not have any other zero; on 
the other hand, when the hyperbolas cross in only two 
points, the Jost determinant has two other zeros, which 
have to form a mutually- conjugated complex pair, as seen 
above. This last case corresponds to a resonance, as il- 
lustrated by Fig. [IJc), where the hyperbolas only have 
two intersection points in the (A, p)-plane and a pair of 



complex roots appears in the complex k and p planes. 
The potential thus has either zero or one resonance. The 
intermediate case of three intersection points for the hy- 
perbolas [Fig.rjjb)] corresponds to the presence of a mul- 
tiple root of Eq. (j32|) , which lies in an unphysical sheet 
(Imfc < 0, Imp > or Irak > 0, Imp < 0) of the Riemann 
energy surface; this case corresponds to a transition be- 
tween a one-resonance and a two- virtual-state situation. 

One sees that the parameters ca\ and a 2 determine the 
position of hyperbola (|34b[) and, hence, the number of 
bound states rt& (0, 1 or 2) and of resonances n r (0 or 
1). Let us now determine, for fixed values of (3 and A, 
the domains in the plane of parameters A = (01,02) 
with constant values of rib and n r . To find domains in A 
where system (|34[) has two complex conjugated roots (one 
resonance), we consider the case where the hyperbolas 
have a common tangent point, as illustrated by Fig.fjjb). 
One can see that the decrease of either a\ or 02 leads to 
the disappearance of the resonance, while the increase of 
either a\ or 02 leads to the appearance of the resonance. 
We define the parametric curves [ai(Ao, po), ck2(Ao, po)] 
in plane A by shifting the tangent point (Xq,pq) along 
the hyperbola p 2 — A 2 = A. These curves limit domains 
in A with either zero or two complex roots. To find them, 
we use the two conditions corresponding to the common 
tangent point (Ao,po) 



Po - 

dp 
dX 



1 



A + a\ 



A=A 



(Aq + OLl)< 



= ±- 



A 



VAf+A 



(35a) 
(35b) 



The upper signs correspond to Ao < (tangent point in 
the second quadrant) while the lower signs correspond to 
Ao > (tangent point in the fourth quadrant). We can 
solve system (|35[) with respect to a\ and o^: 



ai(Ao) = ± 



«2(A ) 




(36a) 



M 7I +sign(A )^A 2 + A.(36b) 



"(Ao' + A)! 

It should be noted that the Schrodinger equation with 
the Cox potential has the following scale invariance: 

ai.2 -> 7"i,2, A^7 2 A, (37a) 
«i,2 -> 7 K i,2, /3^7/3, (37b) 
r —>■ r/7, (37c) 

which leaves A^ = A//3 2 invariant. Hence, we may put 
A = 1 without losing generality. This choice is equivalent 
to measuring energies in units of A. It is convenient to 
express equations ()36j) in terms of dimensionless variables 
ati/0, A d = A//?, Ao -> X /f3: 
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Figure 2: Parametric curves in terms of dimensionless pa- 
rameters defined by Eqs. (|38p in plane A = (5i,a2) for 
A/f3 2 = 1.2. The left-hand-side curves correspond to the 
lower signs in the equations, while the right-hand-side curves 
correspond to the upper signs. The number of resonances, 
n r , is indicated in each domain of the plane. 



These four solutions [taking into account sign(Ao)] can 
be considered as four parametric curves in plane A = 
(ai/P,a2/P), which separate the plane in five regions 
(one inner region and four outer regions, see Fig. [2]). 

In the inner region, the Jost determinant has two com- 
plex roots fei ( 2 = ±k r + iki and, hence, these values of pa- 
rameters ati, oti correspond to one resonance (n r = 1). In 
the four outer regions, the Jost determinant has purely- 
imaginary roots, hence n r = 0. The curves in Fig. [2] tend 
asymptotically to straight lines which are defined as the 
limits for Ao — » and Ao — * ±oo. As a result, one finds for 
all branches two horizontal asymptotes a?/ (3 — ±\ZA^ 
and three slant asymptotes defined by 02/ f3 = — oci/0 
(for the curves in the second and fourth quadrants) and 
CJ2/73 = —oil/ (3 ± 2 (for the curves in the first and third 
quadrants, respectively). 

Consider now the case where the hyperbolas cross at 
the point Ao = 0, po = vA [see the thin dashed lines in 
Fig. [Ha)]. After a small decrease of either ct\ or ct2, the 
number of positive roots, i.e., of bound states, increases 
by one unit. Hence, assuming Ao = and po = VA in 
system ([Ml) , we get the curves 



at (a 2 + VA\ - f3 2 = 0, 



(39) 



which define three domains in the plane of parameters A, 
where Eqs. (|34[) have different number of positive roots 
(see Fig. [3j) . One can directly check that the number rib 
of bound states may be calculated as a function of the 
parameters as 



Figure 3: Curves (|39[) in plane A in terms of dimensionless 
parameters for A//3 2 = 1.2. The number of bound states, rib, 
is indicated in each domain of the plane. 
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-7.5 
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a 2 /P 
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1 



l+o(/i-l)^2, 



(40) 



Figure 4: Regions of the A-plane with different numbers of 
bound states and resonances, (nb,n r ), for the Cox potential 
with A/P 2 = 1.2. 



where the quantities 

Ii = sign (j3 2 - arVA — aia 2 j ■ 1, (41a) 

h = sign(a 2 + VA) ■ 1 (41b) 

may be considered as invariants. For rib = 0, one has 
7i = — 1 and 1% = 1; for n b = 1, one has I\ = 1 and 
I 2 = ±1; for rib = 2, one has I\ = h = — 1- 

Let us now summarize our findings on the number of 
bound states and resonances of the 2x2 Cox potential, 
by combining Figs. [2] and [3] in Fig.[4j where both n b and 
n r are given for all the possible regions of plane A. The 
border lines of these regions, as already discussed, cor- 
respond to the parametric curves defined by Eqs. (|36|) . 
(|3"5|). and to the curves given by Eq. (|3l?|) . From the 
asymptotic behavior of these curves, it is easy to see the 
global structure of the zones. For instance, for the case of 
two bound states, the hyperbolas in Fig. Q] have to have 
four intersection points, which implies that no resonance 
is present. This is the reason why the boundary lines 
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between the zones of bound and resonant states do not 
cross in the lower-half A-plane. Moreover, one can see 
that the topological structure of these zones does not de- 
pend on a particular choice of a parameter A<z = A//3 2 . 
A change of this parameter only leads to a deformation 
of zones, namely, the distance between horizontal asymp- 
totes changes, but does not make any new intersection 
point or new boundary line appear. 

The case of (3 — 0, = oo corresponds to uncoupled 
channels. In this case there are no resonances. Only 
bound or virtual states located in different channels may 
appear (see Sec. HV] A.) 

Up to now, we have excluded the factorization energy 
from our analysis because Eqs. ([34]) are independent of 
ki,2- We will now give a geometrical analysis of condi- 
tions (|25[) , that have to be imposed on parameters k,\ and 
k-2 to warranty a regular character of the Cox potential 
Vr > 0. Let us notice that condition (|25b[) . 

(Ki+ai)(K 2 + a 2 )-/? 2 >0, (42) 

is nothing but Eq. (|34b[) for \ = Ki, p = n%, and a 
shifted value of parameter 1 . Actually, since the left- 
hand side of Eq. (|4"2")l should be positive, there exists a 
positive number C (depending on the set of parameters) 
such that 

(ki + (*i)(k2 + a 2 ) - (f3 2 +C) = 0, C>0. (43) 

This equation represents a hyperbola with the same 
asymptotes as hyperbola (I34bj) but with a larger distance 
between its branches. We thus conclude that the per- 
mitted values of K1.2 are determined by the intersection 
points of hyperbolas (j2~4")l and (|4"3"|) for arbitrary values 
of C, with the additional condition (|25a[) : K\ > —a.\. 
Figure [5] helps to visualize these conditions: the allowed 
values of the factorization energy correspond to the val- 
ues of K1.2 given by the upper-right intersections between 
the solid and the thin-dashed hyperbolas. In Fig. [5J the 
solid and bold-dashed hyperbolas can also be considered 
as representing Eqs. (j34)k the bound-state energies then 
correspond to their (0, 1 or 2) intersections in the first 
quadrant. The allowed values of ki j2 should thus be 
larger than the largest values of A, p corresponding to 
a bound state. The necessary and sufficient condition for 
a regular potential can thus be simply stated as: the fac- 
torization energy should be negative and lower than the 
lowest bound-state energy, if any. 

C. Inversion of the zeros of the Jost-matrix 
determinant 

To solve a realistic two-channel scattering inverse prob- 
lem, it is necessary to express the Cox potential in terms 
of physical data such as the threshold energy, bound-state 
energies, resonance energy and width, or scattering data. 
While the threshold energy explicitly appears in the ex- 
pression of the Cox potential as parameter A, the other 
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Figure 5: Graphical representation of Eqs. (|24|l (solid lines) 
and (|43|) (bold dashed lines for C = 0, thin dashed lines for 
C > 0). The permitted values of ki,3 correspond to the 
upper-right intersections between the solid and thin-dashed 
hyperbolas. After the substitution «i — > A, K2 — * p, the solid 
lines also represent Eq. (|34ap and the bold dashed lines also 
represent Eq. ([34b}, the first-quadrant intersections of which 
correspond to bound states. 



data are directly related to the positions of the zeros of 
the Jost-matrix determinant, as seen above. Ideally, one 
would thus like to directly express parameters ot\, 02, (3, 
and £, which define the Cox potential, in terms of the 
roots of Eq. ([32]) . Certainly, there exist general formulas 
for the roots of the fourth-order algebraic equation (|32p , 
but they are very involved and cannot help much in re- 
alizing the above program. Therefore, we propose here 
an intermediate approach. Two of the roots of Eq. |32|) 
happen to be easily expressed in terms of parameters cti 
and Q 2 - Once two roots are fixed, Eq. (j32|) reduces to a 
second-order algebraic equation for the two other roots, 
thus providing an implicit but rather simple mapping be- 
tween the roots of Eq. ([3"!?]) and the set of parameters. 

Let us denote by (k±,pi) and (^2,^*2) t w0 zeros of 
f(k,p). This imposes some restrictions on the parame- 
ters a\ and a.2- They are not independent anymore from 
the other parameters but should be found as functions 
of fci,2 ? Pi, 2 5 and /? from the system of two equations ob- 
tained from Eq. (|31b|) , written for k = ki, p = pi and for 
k = fc 2 , p — P2 ■ This system reads 

(fci +ia 1 )(p 1 +ia 2 ) +(3 2 = 0, (44a) 
(fc 2 + ia 1 )(p 2 + ia 2 ) + (3 2 = 0. (44b) 

After eliminating parameter 0*2 from system (|44p . one 
finds a second-order algebraic equation for a±, 

a 2 1 -a 1 i(k 1 + k 2 )-k 1 k2+l3 2 ^=0, (45) 

with 

A fe = fca-fci, (46a) 

A p = P2-P1, (46b) 
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from which follows the two possible choices 
i(ki + k 2 ) ± 



0(2 



A 2 - 4/?2A fc /A p 



,(47a) 



i(pi+P2)T\/-A2-4/32A p /A fe .(47b) 



The upper (resp., lower) sign in Eq. (|47b[) corresponds 
to the upper (resp., lower) sign in Eq. (|47a|) . The values 
of k\ t 2 and pi j2 should be chosen so as to warranty the 
reality of parameters ax 2- 



Now, two roots ki and k2 of the fourth-order alge- 
braic equation (|32|) are fixed. Therefore, the two re- 
maining roots &3 and ki are solutions of a second-order 
equation Q 2 (fc) = 0, where polynomial Q2{k) is the 
ratio of the polynomial appearing in Eq. (|32[) and of 
P2(k) = k 2 — k(k2 + ki) + fefei, i.e., explicitly 



fc 4 + iaxk 3 + a 2 k 2 + ia 3 k + a 4 = P 2 (k)Q 2 {k) . 
From here we find 



Q 2 {k) = (k + ioixf + k{k 2 + fa) + (2iai + k 2 + k 1 ){k 2 + fa) + a\ - A - fak 2 



(48) 



r 



and, hence, 



k'4 



TU-A 2 -A(3 2 A k /A p + ^D k 



TiJ-A 2 - A(3 2 A k /A p - ^D k 



(49a) 
(49b) 



where D k = Af +4/? 2 ^+4fc 2 fci. The sign before the first 
square root in Eqs. (|49|) should be chosen in accordance 
with the sign in Eqs. (|47|) for parameters ai. 2 . To find 
P3,4, we do not need to solve any equation. We simply 
notice that the equation det-F(fe) = is invariant under 
the transformation k <-* p, a\ «-> a 2 , A <-» —A. This 
means that, being transformed according to these rules, 
Eqs. (j49|) give us the p values: 



1 

P3= 2 



Pi= 2 



Ti^-A 2 p - Ap 2 A p /A k - ^D p 
TiJ-A 2 p - 4/? 2 A p /A fc + y/Dl 



(50a) 
(50b) 



where D r , 



, p -A 2 +A(3 2 ^+4p 2 p 1 . 
Let us now consider several important examples. Ac- 
cording to the analysis of Subsec. IIIIB1 only two essen- 
tially different possibilities exist, namely, a system de- 
scribed by the Cox potential may have either one reso- 



nance or no resonance. 



1. One resonance 



Assume Eqs. (|3 1 1) have two complex roots and let us 
define their first-channel components as 



ki 

k 2 



hip ~\~ ik'i i 



(51a) 
(51b) 



where k r and ki are real. Let us assume the real part 
k r to be positive to fix ideas. The imaginary part ki, on 



the other hand, can be either positive or negative. Let us 
write the corresponding energies, k\ 2 , as E r ±iEi, where 
we also assume E 1 , positive to fix ideas (which means that 
the upper sign corresponds to k\ or fc 2 , depending on the 
sign of ki). We would like to choose as parameters the 
threshold difference A, as well as the real and imaginary 
parts of the resonance complex energy, E r , Ei. As seen 
below, these can correspond to physical parameters of a 
visible resonance in some (but not all) cases. In terms of 
these parameters, k r and ki are expressed as 



k * 
r V2 



El + E 2 - E r 



h = ± 



1 



V2 



-1/2 



1/2 



El + E 2 - E r 



(52a) 
(52b) 



In the second channel, the roots corresponding to /ci ]2 
can be found from the threshold condition l|28[). They 
are given by 



Pi 

P2 



Vr + iPi, 
-Pr + IPi 



(53a) 
(53b) 



with 



Pr = 



1 

V2 
Ei_ 
V2 



(E r - A) 2 +E r — A 



1/2 



,(54a) 



E 2 + (E r - A) 2 +E r — A 



-1/2 



.(54b) 



The upper (resp., lower) sign in Eq. ()52b[) corresponds to 
the upper (resp., lower) sign in Eq. (|54bj) . which means 
that, for a given zero, the signs of ki and pi are opposite. 
This can be understood from the first column of Fig. [1] 
a resonance appears when the two hyperbolas are tan- 
gent to each other, which can only happen in the second 
and fourth quadrant, where A and p have opposite signs. 
Moreover, Eqs. (|52a[) and (|54aj) show that, for a given 
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zero, the signs of k r and p r are also opposite. This im- 
plies that, for the Cox potential, the complex resonance 
zeros (or scattering-matrix poles) are always in opposite 
quadrants in the complex k and p planes, as illustrated 
for instance by the complex zeros in Fig. HJc). 

This has important consequences for physical applica- 
tions: for a resonance to be visible, one of the correspond- 
ing zero has to lie close to the physical positive-energy 
region, i.e., close to the real positive k axis and close to 
the region made of the real positive p axis and of the 
positive imaginary p interval: [0, iyA]. Consequently, 
the only possibility for a visible resonance with the Cox 
potential is that of a Feshbach resonance, only visible in 
the channel with lowest threshold, with an energy lying 
below the first threshold. At higher resonance energies, 
the corresponding zero is either close to the ft-plane phys- 
ical region (and far from the p-plane one) or close to the 



p-plane physical region (and far from the fc-plane one); 
it cannot be close to both physical regions at the same 
time, hence it cannot have a visible impact on the cou- 
pled scattering matrix. Here, we limit ourselves to the 
case of a visible resonance, which is the most interesting 
from the physical point of view. It corresponds to the 
lower signs in Eqs. (|52bj) and (I54bj) . with a resonance en- 
ergy E r such that < E r < A, and a resonance width 
T = 2Ei such that E { < E r . 

To get a potential without bound state, we choose the 
upper signs in Eqs. (|4"T)) , which leads to 

ai = + k r [-1 - p 2 /(k rPr )] 1/2 , (55a) 

1/2 

a 2 = -Pi - p r [-1 - (3 2 /(k r p r )] . (55b) 
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Figure 6: The Cox potential with one visible resonance of energy E r = 0.4 and width F — 0.02, for A = 1 and = 0.1 (first 
row, solid lines for fin and V22 + A, dashed line for V12), with the corresponding partial cross section (second row) and phase 
shifts (third row) for (a) ki = 0.5; (b) «j = 0.7; (c) ki = I. 



From here, we see that, for non-zero values of the pa- 
rameters k r and p r (which have opposite signs), the cou- 
pling parameter (3 cannot be infinitesimal: because a\ 
and Q2 have to be real, (3 is restricted to satisfy the in- 



equality (3 > \J—k r p r . Now, using Eqs. (|49|) and (f50|) it 
is easy to find the two remaining roots 



k 3A = —iy/—kl - (3 2 k r /p r ± iJ ft? - (3 2 p r /k r , (56a) 
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P3,4 = -i^/-Pr ~ P 2 Pr/k r T iy P 2 ~ P 2 k r /p r ■ (56b) 

To get a potential with one bound state at energy — A^, 
we choose the lower signs in Eqs. (j47|) . We then get 
for ks,{(3) an expression similar to Eq. (I56a|) . from which 
the value of (3 can be found by solving the bi-squared 
equation 

k 3 (P)=i\ b . (57) 

Let us now choose explicit parameters. First, without 
loosing generality we may put A = 1 (see Subsec. UlI Bl) . 
To get a visible resonance, we put E r = 0.4, = 0.01 
(which corresponds to a resonance width T = 0.02), and 
= 0.1. Using Eqs. <P J) and ([SB ) , one finds a x = 0.76938 
and «2 = —0.766853. The factorization energy, £, is not 
constrained in this case: it just has to be negative. The 
Cox potential with one resonance and no bound state for 
different values of n\ is shown in the first row of Fig. [6] 
The diagonal elements of the potentials, V\x and V22 + A, 
are plotted with solid lines, while V\% is plotted with 
dashed lines. Parameter ki is responsible for changing 
the range of the potential, as shown by Eqs. (|26|) . and 
hence the scattering length (see discussion below). The 
second row of this figure shows the corresponding partial 
cross sections, where the resonance behavior is clearly 
seen, as well as the evolution of the low-energy cross sec- 
tion, which is related to the scattering length. The last 
row of Fig. [5] shows the corresponding phase shifts for 
the open channel, where a typical Breit-Wigner behavior 
(see e.g. Ref. [6j) is seen for the resonance, as well as the 
evolution of the zero-energy phase-shift slope, which is 
also related to the scattering length. 




Figure 7: The Cox potential (solid lines for Vu and -322 + A, 
dashed line for C12) with two bound states at energies E\ = 
-0.01 and E 2 = -2.25, for A = 1, /? = 0.1 and Ki = 1.51. 
The left (resp., right) graphic corresponds to the upper (resp., 
lower) signs in Eqs. (I47|l . 
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Figure 8: Iso-energy curves defined by Eq. (|58[) (dashed lines) 
in the plane of parameters (01,02) for Eb = —4, A = 1, 
f3 = 2, together with the curves defined by Eq. (|39p (solid 
lines) . 



2. Two bound states 

Let us now construct a Cox potential with two bound 
states, and hence no resonance (see Fig. 0}. We choose 
k\ = O.li and A)2 = 1.5* for these bound states and, as in 
the previous example, we put A = 1 and (3 = 0.1. Wc 
thus havepi = yLOlz andj>2 = V3.25i, which defines A p 
in Eqs. (|4T[) . Choosing the upper signs in these equations, 
we find a± = —0.112649 and 012 = —1.79557, while for 
the lower signs, we get ol\ = —1.48735 and 012 = —1.0122. 
The corresponding Cox potentials are shown in Fig. [7] 



3. One bound state 



Figure [5] gives a graphical representation of Eq. (f5~5|) in 
the plane of parameters (ai,a2): the dashed curves cor- 
respond to all the Cox potentials that have a bound state 
at the same energy Ej, = —4. Figure [5] also shows that 
these iso-energy curves ([55| can be obtained by a shifting 
of the curves separating regions of plane A with different 
number of bound states, defined by Eq. (|39|) (solid lines). 

Let us finally summarize these possible inverse prob- 
lems in Table [U The free parameters allow either for 
isospectral deformations of the potential or for fits of ad- 
ditional experimental data as, e.g., scattering lengths. 
This possibility will be used in the next section on atom- 
atom systems. 



If only one bound state E b = — is known, we may 
consider the position of the second zero as a free param- 
eter. Besides, one can solve system (|34j) at A = A& with 
respect to «2, which leads to 

a 2 = y/A + Ag. (58) 



D. Low-energy scattering matrix 

In this section, we analyze the S'-matrix given by 
Eq. ((3]) for energies close to the lowest threshold, the 
energy of which we have chosen equal to zero. From 
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Table I: Possible mappings between some experimental data 
and the Cox potential parameters. 



Experimental 


Fixed 


Free 


Restrictions 


data 


parameters 


parameters 




A , E r , Ei 




Kl,(3 


/3 > ^/ — k r p r 


A , E b = — A^ ,E r ,Ei 


«1 ,Q!2,/3 


Kl 


Kl > Af, 


A , -El,2 = — A? i2 


ai , a 2 


Kl,/? 


ki > A 2 > Ai 


A , E b = -Ag 


a 2 




Kl > X b 



Eqs. (f2"9"| and (f3T)]) . one finds the Cox-potential S-matrix 
1 



S(k,p) = 



f(k,p) 



~2if3^/kp 



m -p) 



(59) 



When the second channel is closed, i.e., for energies < 
E < A, the physical scattering matrix is just a function 
S{k,p), which coincides with the first diagonal element 
of S'-matrix (l59l) . It reads 



S(k,p) 



k + iK\ [i(k - iai)(VA - k 2 + a 2 ) - 1 
k - iKi [i{k + zai)(VA - k 2 + a 2 ) + /3 2 



(60) 

From here one finds the scattering amplitude 
A(k) = [S(k) - l]/2ik, which reads 



A(k) = - 



h + VA-fc 2 )(ai-«i)-/3 2 



i (k - iKi) [i(k + iai) (VA - k 2 + a 2 ) + (3 2 
and the scattering length a = —A(0), which reads 



(61) 



1 



Oi 2 



(62) 



From the argument of S(k) — e 2lS ( k \ one deduces the 
phase shift S(k), which reads 



k k(VA~k^ 

o(k) = arctan h arctan — 

hi (3 2 



0-2) 



ai (\/A - k 2 + a 2 ) ' 
(63) 

One can check on Eqs. (|6"2"]) and ([6^)) that the scattering 
length is the slope of the phase shift at zero energy, as it 
should be. Note that Eq. ([63]) is equivalent to 



k cot S(k) 



ap(k)Ki 



ki - ap 



(64) 



where ap(k) = ai — [3 2 / (VA — k 2 + a 2 ) . In the un- 
coupled case (/3 = 0), this expression reduces to the 
phase shifts of the simplest Bargmann potential (see 
e.g. Ref. [H|]), which depends on the parameters ki and 
cib = a/3=o = a i- Therefore, the Cox potential may 
be considered as a coupled-channel deformation of the 
Bargmann potential, resulting in an energy dependence 
of one of its parameters, ag. 



The scattering length is an important physical quan- 
tity. In many-body theories for instance, it is often used 
to describe interactions in the s-wave regime. Let us thus 
study into detail the scattering length of the Cox poten- 
tial, as given by Eq. (|62[) . When considered as a function 
of ai j2 , it has a singularity located at the boundary of 
the single-bound-state region provided by Eq. ([39]) . Such 
infinite values of the scattering length happen when a 
zero of the Jost determinant, which corresponds to an S'- 
matrix pole, crosses the first threshold: a bound state is 
then transformed into a virtual state, in agreement with 
the general theory [l2|] ■ We can analyze the sign of the 
scattering length for different numbers of bound states rib 
given by Eq. ([40]) . by considering the indices Ii >2 given by 
Eqs. flT]) . First, we remark that Ii = —1 for n b = 2,0, 
and Ii — 1 for rib = 1- Then, since 1% = — 1 for rib = 2, 
the scattering length is positive when the Cox potential 
has two bound states. Next, for n b — we have I 2 = I; 
the two contributions of the scattering length then have 
different signs. As a result, at fixed ai, a 2 , (3, and A, 
one can get both positive and negative scattering lengths 
by varying only ki. Similarly, when I 2 = I (one-bound 
state), the scattering length may only be positive. In 
contrast, for I2 = —1 it may be both positive and nega- 
tive but for fixed values of ai, a 2 , (3, and A, it becomes 
negative for large enough k%. 



IV. TWO-CHANNEL MODEL OF 
ALKALI-METAL ATOM-ATOM COLLISIONS IN 
THE PRESENCE OF A MAGNETIC FIELD 

A. Magnetic Feshbach resonance 

Ultra-cold collisions of alkali-metal atoms play a key 
role in applications of laser cooling such as Bose-Einstein 
condensation and BEC-BCS crossover. The analysis of 
such experiments is commonly based on the coupled- 
channel method [l8j . i.e., on solving numerically a set 
of coupled differential equations. 

In this paper, we reduce the low-energy scattering 
problem of two alkali-metal atoms to an effective two- 
channel problem with a single Feshbach resonance, as in 
Ref. 0. The model consists of a single closed channel 
Q containing a bound state, which interacts with the 
scattering continuum in the open channel P, so that the 
whole scattering problem is reduced to the two-channel 
scattering described by the 2x2 Hamiltonian 



H = - 



dr 2 



V P {r) V lnt {r) 
V int {r) Vg(r) 



(65) 



where Vp is the uncoupled open-channel potential, Vq 
is the uncoupled closed-channel potential, and potential 
Vint describes the coupling between the open and closed 
channels P and Q. These channels describe atoms placed 
in a magnetic field and occupying different energy sub- 
levels which can be shifted with respect to each other 
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with the change of the magnetic field (Zeeman effect). 
For each value of the magnetic field, the zero of energy is 
chosen as the energy of the dissociated atoms in channel 
P. 

Even in the simplest case of a homogeneous magnetic 
field, the potential-energy matrix of Hamiltonian (|65|) de- 
pends on the magnetic field. We will assume that the 
external field changes slowly enough so that we can take 
advantage of the adiabatic approximation, assuming that 
the stationary Schrodinger equation may be applied for 
describing the scattering process and the magnetic field 
enters the Hamiltonian as a parameter only. Moreover, 
the known observation that, when the scattering length 
is much larger than the range of the interaction, the gen- 
eral behavior of the system is nearly independent of the 
exact form of the potential [19j, suggests us to use the 
Cox potential with large scattering length for describing 
the interatomic scattering. We thus replace the potential 
matrix in Eq. (|65p by the Cox potential. In this case, the 
parameters of the Cox potential should carry a depen- 
dence on the magnetic field. Below, we show that, to 
get a good agreement with available experimental data, 
it is sufficient to impose a linear field dependence on the 
threshold difference A only, keeping all other parameters 
field independent. Thus, inverting known scattering ex- 
perimental data, one can find all the parameters defining 
the Cox potential, obtaining in this way a simple ana- 
lytical model of the atom-atom scattering process in the 
presence of a magnetic field. 

The position of the highest bound (or virtual) state is 
crucial in describing the resonance phenomena of inter- 
atomic collisions. In an s-wave single- channel system, the 
scattering process becomes resonant at low energy when 
a bound state or virtual state is located near the thresh- 
old, a phenomenon known as "potential resonance" . In a 
multichannel system, the incoming channel (which is al- 
ways open) may be coupled during the collision process 
to other open or closed channels, corresponding to differ- 
ent spin configurations. When a bound state in a closed 
channel lies near the collision energy continuum, a Fesh- 
bach resonance P, [|J may occur, giving rise to scattering 
properties that are tunable by an external magnetic field. 
In Ref. 0] , some interesting examples of the interplay be- 
tween a potential resonance and a Feshbach resonance are 
considered. Below, we adjust the analytically-solvable 
model based on the Cox potential for describing the same 
phenomena. 

Typically, the coupling between the closed and open 
channels is rather small; we thus consider first an uncou- 
pled limit of the Cox potential, i.e., Vi n t{r) — > 0, which 
corresponds to j3 — > 0. In this case, the Jost determi- 
nant (|30[) has the following zeros: 



(66) 



and 



The energies of these unperturbed (i.e., with zero cou- 
pling) states (called bare molecular states in Ref. 0]) are 



Ex 



and 



E 2 = -ai 



(68) 



(69) 



It should be noted that in this case E\ belongs to channel 
P while E2 belongs to channel Q. Hence, ai is associated 
with the potential resonance, while 012 is associated with 
the Feshbach resonance. Due to the Zeeman effect, the 
difference between the thresholds is a linear function of 
the magnetic field, 



MB) 



Hmag{.B — Bq), 



(70) 



where Bq can be arbitrarily chosen in the domain of inter- 
est and Ao is the value of the threshold corresponding to 
Bq- If «i,2 < and the coupling is absent, then the two 
bound states cross at A = a\ — a\. Note that E2 crosses 
the threshold at A = a\. When there is a coupling be- 
tween channels, the levels E\ and E2 avoid crossing (see 
below) . 

Let us consider the behavior of the scattering length in 
the presence of the Feshbach resonance. It is described 
by the following formula 



1 



r, 



B-Bo 



(71) 



Here, Bq is the position of the magnetic Feshbach reso- 
nance and is its width (in terms of magnetic field) . 

In particular, Eq. (j6"2")) shows that such an infinite value 
of the scattering length occurs for the Cox potential at a 
threshold Aq defined by: 



/3 2 — ai« 2 

Oil 



(72) 



Let us now assume for the Cox potential a threshold dif- 
ference given by Eq. (|T0[) with such a value of Ao. Ex- 
panding Eq. (f6"2"|) near this resonance, it is easy to get 
Eq. (|7ip and a simple expression for the width of the 
resonance for the Cox potential. 



OL\ — K\ 



OL\K\ 



(73) 



2 [1 + o (A - A )] ki ( \[M + a 2 ) 



(ai-«i)(Ao-A) 



2KiV^o"( v /An" + a 2 ) 

[i-rnag [ol\ — K\) 



(74) 



P2 



-ia%. 



As shown in Ref. Q , the background scattering length 
(67) ab g is due to the open-channel potential. When there is 
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a bound state or virtual state close to threshold, it can 
be further decomposed as a sum of two contributions: 
a standard potential part, which depends on the poten- 
tial range, and a potential-resonance part, which depends 
on the bound/virtual-state energy. This decompostion 
clearly appears in our model: the background scattering 
length corresponds to a large magnetic field B in Eq. (PTl"]) 



above B = 0.112, when the imaginary part of the zero in 



or to a large A in Eq. ([52")) , namely, at, g 
yields: 



lim a, which 



A^oo 



a bg 



1 

«1 



1 



(75) 



The same result comes from Eq. ([To]) . It should be noted 
that a large threshold difference A — > oo effectively cor- 
responds to a small coupling (3 — > 0. From here we find 
another relation ab a — lim a which leads to the same 

13^0 

background scattering length. 

In this formula, the first term is proportional to 
the parameter which defines the range of the open- 
channel potential [see Eqs. ([26]) ]: it may thus be con- 
sidered as the standard potential part of the background 
scattering length. The second term is associated with 
the P-channcl bound (or virtual) state in the uncoupled 
limit. Hence, it may be interpreted as the potential- 
resonance part of the background scattering length. Let 
us further consider two different possibilities giving rise 
to a large (either positive or negative) background scat- 
tering length. 



Interplay between a bound state and the 
Feshbach resonance 



The first possibility occurs when the highest bound 
state is located near the threshold, i.e., when a,\ < 0. 
In Fig. [9j we show energies as functions of the magnetic 
field when channel P has a bound state just below the 
threshold, for 





Ol 



0.05, 

-A 6 = -0.103, 

-0.5, 

1. 



(76a) 
(76b) 
(76c) 
(76d) 



We are using arbitrary units and choose A(B) = 
0.35 - B in Eq. ([TDD. The bare bound states [fl = 0) of 
the P and Q channels are indicated by the dashed hor- 
izontal and slanted lines respectively. These bare states 
are Hamiltonian eigenstates in the uncoupled P and Q 
subspaces. The bare Q-channel bound state crosses the 
P-channel threshold (solid horizontal line) at B = 0.1. 
The dressed states are represented by solid lines and dis- 
play an avoided-crossing behavior For low fields, the 
model has one bound state and one Feshbach resonance, 
the energies of which are close to the bare-state energies. 
The energy of the Feshbach resonance becomes negative 
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Figure 9: Energies of bare (dashed lines) and dressed (solid 
lines) states as functions of the magnetic field B for the Cox 
potential defined by parameters (|76[) . The transition between 
a Feshbach resonance, a virtual state, and a bound state is 
shown in the inset. 




Figure 10: The Cox potential defined by parameters (|76[) for 
B — 0.1; Vp and Vq + A are represented by solid lines, Vint 
by a dashed line. 



the complex k\ plane becomes larger than its real part. 
At B = 0.12, these complex zeros collapse and transform 
into two virtual states (purely imaginary zeros), which 
corresponds to the discontinuous slope in Fig. [5] With 
increasing magnetic field, one of these virtual states (rep- 
resented in Fig. [9]) gets closer to threshold, while the 
other one (not represented in Fig. [5] as it does not af- 
fect the low-energy scattering properties) goes away. At 
Bo = 0.124, the virtual state crosses the threshold and 
becomes a bound state; the scattering length thus goes 
through infinite values at that field: this is the magnetic- 
Feshbach-resonance phenomenon itself. 
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Above Bq, the model has two bound states, the energies 
of which tend to the bare-state energies when the field 
increases. 

Following Ref. 0, we stress that, although the behav- 
ior of the dressed states shows some resemblance with the 
two-level Landau-Zener description, this model does not 
include the threshold effects shown in Fig. [S] and, hence, 
cannot be used to properly describe the interplay between 
a potential resonance and a Feshbach resonance. With 
respect to Ref. Q , our model displays a slightly more so- 
phisticated behavior for the state energies (compare our 
Fig. [9] with their Fig. 4). A more significant novelty of 
our description is the direct knowledge of the coupled- 
channel potential corresponding to these energies. This 
potential is shown in Fig. [TU] for B = 0.1. The potential 
form factor changes slowly with the change of the mag- 
netic field, which is mainly responsible for the variation 
of A. 

The value of k\ chosen in Eq. (|76d[) is arbitrary. How- 
ever, the necessary and sufficient condition to get a 
Cox potential without singularity imposes then that the 
bound-state energies of the model should be larger than 
— 1. Figure [5] shows that this condition will be satisfied 
for a limited range of magnetic field only. For higher 
fields, a larger n\ should be chosen. The phase shifts of 
the same Cox potential, as well as a graphical representa- 
tion of Eqs. ([34]) , are shown in Fig. [TT] for different values 
of B. The first and the last columns correspond to a large 
positive background scattering length (ab g ~ 1/A& « 10), 
due to a bound state close to the threshold. Physically, 
this occurs for the 133 Cs atom-atom interaction |20|, for 
instance. FigureQTJb) illustrates the case where the scat- 
tering length is close to zero. The calculation or measure- 
ment of the zero of the scattering length plays an impor- 
tant role in determining the resonance width [2l| . The 
phase-shift behavior for the virtual state and bound state 
close to threshold is shown in Figs.QTJc) and [TlTd). re- 
spectively. In this case, the scattering length is very large 
and its sign changes while the energy of the zero of the 
Jost-matrix determinant crosses the threshold. Recalling 
that the intersection points in the graphical representa- 



tion of Eqs. ([34]) , shown in the second row of Fig. QTJ give 
the positions of bound and virtual states, one may estab- 
lish a correspondence between the second row of Fig. [TT1 
and the motion of the corresponding zeros in the complex 
plane described above. 



Interplay between a virtual state and the 
Feshbach resonance 



Another interesting possibility occurs when there is a 
virtual state close to the threshold, i.e., when ot\ > 0. 
This is the case of the 85 Rb atom-atom interaction, for 
example. We will use rubidium scattering data 0, H2| in 
this example, and work with units h = 2/i = 1, where \x 
is the reduced mass of the two atoms. The length unit is 
chosen as the Bohr radius ao; energies arc thus expressed 
in units of a^ 2 . According to Ref. 0], the bare virtual 
state is located at A„ = —1.78 • 10~ 3 aQ , but this value is 
associated with the model they used in their calculations. 
We just consider X v ~ — 10 _3 a ( ^ 1 and set Eq. (175)) as a 
constraint between a.\ = —X v and K\. In order to fit 
the scattering-length behavior (fTTj) with ab g = — 443 ao, 
B = 15.5041 mT and r B = 1.071 mT, we use Eq. ((62]) . 

The value of [3 defines, in particular, the position of 
the Feshbach resonance, i.e., the magnetic field Bq for 
which the bound state crosses the threshold. According 
to Eq. (|39|) , one has 



a>i a 2 



(77) 



where Ao is the value of the threshold corresponding to 
Bo. The value of defining the width of the Feshbach 
resonance Tb, should be found from the condition a(B + 
T B ) = 0. Then, according to Eq. (|6"2"j) . we find 



ax 



a 2 



^A(B + T B ) 



Ki 



(78) 
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Figure 12: The Cox potential describing the Feshbach resonance in 85 Rb, defined by parameters (|79[1 . plotted at £> = 14.5 
mT (A = 0.0590363 a^ 2 ). 
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Figure 13: Solid line: Feshbach-resonance scattering 
length (I71|) . Dots: Cox-potential scattering length (|62|l for 
the parameters (|79l) . 



where Ao = 2471.386 MHz and ^ mag = 
-36.4 MHz/mT 0. To get that value of A , we 
have used the known value of the threshold at zero 
magnetic field [22| and assumed that Eq. ([70)1 is valid 
down to that field. 

From Eq. (|75)l . we may fix k\ — ai/(l+af> s Ki) at a,b g = 
—443 ao and find the values of all parameters defining the 
potential at the given position of the Feshbach resonance 
and with the given value of the background scattering 



length: 

(3 = 0.0202366 Oq 1 , (79a) 

ax = -X v = 2.2 • 10~ 3 a \ (79b) 

a 2 = -0.239343 a^ 1 , (79c) 

k-l = 0.0866 Uq 1 , (79d) 



k 2 = yjK\ + A = V0.0789668- 0.856899Sa " 1 .(79e) 

In Fig. [13l we show that, with these parameters, 
the Cox-potential scattering length (|62|) reproduces the 
Feshbach-resonance scattering length ([TTj) with good pre- 
cision. The value a\ = 2.2 • 10 _3 a ~ 1 was chosen to get 
a smooth potential Vp without repulsive core. This po- 
tential is shown in Fig. [T^] and, once again, has a form 



factor rather independent of the field, except for the 
threshold. In Fig. 1151 we show the corresponding en- 
ergies as functions of the magnetic field. The bare bound 
state of channel Q is represented by the slanted dashed 
line. The bare virtual state of channel P, which is lo- 
cated at A„ — —2.2 • 10 _3 a " 1 , is not shown in Fig. \TE\ 
The dressed states are indicated by solid lines. When 
B < Bq = 15.5041 mT, there exist both a virtual state 
and a Feshbach resonance, the energies of which tend to 
the bare-state energies for small B. The virtual state 
becomes a bound state at B = Bq (see inset). With 
increasing B, the real part of the resonance energy de- 
creases and at B = 16.657 mT it crosses the threshold. 
Finally, at B = 16.9 mT, the two resonance poles collapse 
and produce two virtual states, one of which stabilizes at 
A„ = —2.2 • 10 _3 a " 1 ( the other one has a much larger 
negative energy and is not represented in Fig. 1151 as it 
does not affect the low-energy scattering properties). 

The behavior of the curves in Fig. [15] is very similar to 
those of Fig. [0J in particular regarding the transformation 
of the Feshbach resonance into a virtual state. The only 
difference between the present case (avoided crossing be- 
tween a virtual state and a Feshbach resonance) and the 
previous case (avoided crossing between a bound state 
and a Feshbach resonance) is that here a virtual state 
transforms into a bound state before the crossing, while 
there a virtual state transforms into a bound state after 
the crossing. Another interesting comparison is between 
our Fig. M and Fig. 5 of Ref. 0; it would be instructive 
to perform a detailed comparison of the two models to 
explain the differences between these two figures. 

As for the interplay with a bound state, Fig. [15] also 
shows some limit on the range of magnetic field on which 
our model can be used: since k,\ is fixed in Eq. (|79d|) 
and the bound-state energy should be larger than — n\ « 
— 0.0075aQ 2 (otherwise the potential becomes singular for 
some value of r), the field should be lower than 24.5 mT. 

The behavior of the phase shifts in the region with 
the resonant and virtual states is shown in the first row 
of Fig. [T3] A similar discussion to that of Fig. [11] can 
be made here, except that here the large positive back- 
ground scattering length results in a large negative slope 
for the phase shift at the origin. 
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Figure 14: Phase shifts and graphical representation of Eqs. (|34p for the Cox potential defined by parameters (I79|l . 
The columns correspond to different values of the magnetic field: (a) B = 14.454 mT; (b) 15.504 mT; (c) 15.854 mT ; 
(d) 19.0 mT. 
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Figure 15: B-dependence of the energies of the bare (dashed 
lines) and dressed (solid lines) states for the Cox potential 
defined by parameters (|79[1 . 



Exactly at Bq — 15.5041 mT, when the bound state 
transforms into a virtual state, the phase shift starts from 
7r/2. The second row of Fig. [TH shows the correspond- 
ing behaviour of the bound- and virtual-state zeros on 
the wave- number imaginary axes, confirming the above 
analysis. 



V. CONCLUSION 

In this work, we have derived the exactly-solvable N- 
channel Cox potential from a supersymmetric transfor- 
mation of the vanishing potential and we have established 
different parameterizations of this potential, as well as a 



necessary and sufficient condition for its regularity. In 
the N = 2 full analysis of the corresponding Jost 

matrix has been carried out. The structure of the zeros 
of the Jost determinant has been presented geometrically 
and a method for controlling the position of the zeros 
of this Jost determinant has been proposed. This has 
led to several examples of Cox potentials with different 
number of bound states and resonance, solving schematic 
coupled-channel inverse problems. 

With ultracold gases in mind, we have also studied the 
low energy S'-matrix and the scattering length of the Cox 
potential. Using independence of scattering properties 
from interaction details in the regime with a large scat- 
tering length, a model of alkali-metal atom-atom scat- 
tering has been constructed. This provides interesting 
exactly-solvable schematic models for the interplay of a 
magnetically-induced Feshbach resonance with a bound 
state or a virtual state close to threshold. 

We consider the development of supersymmetric trans- 
formations as a very promising tool for the multi-channel 
inverse scattering problem and for the construction of 
more advanced exactly-solvable coupled-channel mod- 
els. In particular, iterations or chains of transformations 
might lead to more complicated Jost functions, with arbi- 
trary number of bound states and resonances, hopefully 
still with a tractable connection between potential pa- 
rameters and physical observables. 

As far as physical applications are concerned, atom- 
atom interactions are both very interesting today, due 
to the active research field of ultracold gases, and rather 
simple with respect to supersymmetric quantum mechan- 
ics, as only s- waves have to be considered and as the 
interaction is short ranged (no Coulomb term) . We fore- 
see to apply the present model to other systems present- 
ing these simple features, namely coupled s-wave baryon- 
baryon interactions, with at least one neutral baryon. In 
the longer term, we hope to generalize our method to 
higher partial waves and to Coulomb interactions. This 
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should allow us to construct useful models in the con- 
text of low-energy nuclear reactions, the field which first 
motivated the work of Feshbach P, Q on coupled-channel 
resonances, leading to possible applications in nuclear as- 
trophysics and exotic-nuclei low-energy reactions. 
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